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Abstract: We propose a simple but practical methodology for the quantification of correlation risk in 
the context of credit derivatives pricing and credit valuation adjustment (CVA), where the correlation 
between rates and credit is often uncertain or unmodelled. We take the rates model to be Hull-White 
(normal) and the credit model to be Black-Karasinski (lognormal). We summarise recent work 
furnishing highly accurate analytic pricing formulae for credit default swaps (CDS) including with 
defaultable Libor flows, extending this to the situation where they are capped and/or floored. We also 
consider the pricing of contingent CDS with an interest rate swap underlying. We derive therefrom 
explicit expressions showing how the dependence of model prices on the uncertain parameter(s) can 
be captured in analytic formulae that are readily amenable to computation without recourse to Monte 
Carlo or lattice-based computation. In so doing, we crucially take into account the impact on model 
calibration of the uncertain (or unmodelled) parameters. 

Keywords: perturbation expansion; Green's function; model risk; model uncertainty; credit 
derivatives; CVA; correlation risk 


1. Introduction 

1.1. Model Risk Management 

Much effort is currently being invested into managing the risk faced by financial institutions 
as a consequence of model uncertainty. One strand to this effort is an increased level of regulatory 
scrutiny of the performance of the model validation function, both in terms of ensuring that adequate 
testing is performed for all models used for pricing and risk management purposes and of enforcing 
a governance policy that only models so tested are so used. As is stated in the Supervision and 
Regulation Letter of US Federal Reserve (2011): 

An integral part of model development is testing, in which the various components of a model and its 
overall functioning are evaluated to show the model is performing as intended; to demonstrate that it 
is accurate, robust, and stable; and to evaluate its limitations and assumptions. 

Another concern is model risk monitoring and management. Here, the idea is that, having validated 
models and examined the associated uncertainty, the risk department should monitor and report on the 
risk faced by a financial institution, ideally so that senior management can, based on "risk appetite", 
make informed decisions about model usage policy. According to US Federal Reserve (2011): 

Validation activities should continue on an ongoing basis after a model goes into use to track known 
model limitations and to identify any new ones. Validation is an important check during periods 
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of benign economic and financial conditions, when estimates of risk and potential loss can become 
overly optimistic and the data at hand may not fully reflect more stressed conditions...Generally, senior 
management should ensure that appropriate mitigating steps are taken in light of identified model 
limitations, which can include adjustments to model output, restrictions on model use, reliance on 
other models or approaches, or other compensating controls. 

Here, the notion of best practice is less well established, in particular because different institutions 
adopt different approaches to measuring and reporting model risk. In what probably remains the most 
definitive book 1 on the subject, Morini (2011) asserts: 

You will see that not even among quants there is consensus about what model risk is. 

Indeed, there are currently regular industry events held at which practitioners and managers from 
financial institutions share and discuss their views on current best practice and how this should evolve. 
It is not therefore possible to enforce specific regulatory standards in this area, although regulators do 
take an interest in how banks perform the model risk governance function. 

Central to the task of monitoring and managing model risk or uncertainty is the challenge of how 
to measure it. Current practice tends to be a mix of qualitative and quantitative metrics. While the 
former are easier to implement, the latter are preferable in terms of the level of control that can be 
exercised, particularly if the model risk can be quantified in monetary terms. However, the fact that no 
commonly agreed methodology has emerged and such methodologies as have been proposed tend not 
to lend themselves to implementation by practitioners means that it is not easy to make progress in 
this area. 

A common approach taken by financial institutions has been to consider the reserves taken by the 
finance function to account for model parameter and/or calibration uncertainty as a proxy measure of 
model risk. This is arguably less than satisfactory for a number of reasons, not least that the purpose of 
reserves is to provide a protective buffer against, rather than a precise measure of, model risk. 

The present paper represents a proposed compromise between rigour and practicality to furnish 
model risk metrics against which risk appetite can be compared. We consider specifically rates-credit 
correlation risk in relation to credit derivatives pricing, but it is suggested the methodology may be 
applicable more widely to other types of model risk and a wider class of financial instruments. 

1.2. Layout of the Paper 

We begin in Section 2 by reviewing previous methodologies that have been proposed for the 
quantification of model risk, before formally outlining our own proposed approach. We go on in 
Section 3 to describe the model we shall consider for pricing credit derivatives in the context of 
stochastic interest rates and credit intensity, indicating how, following Turfus (2017a), a Green's 
function solution can be obtained as a perturbation expansion, under the assumption that both of these 
rates are small. In Section 4, the key results of Turfus (2017a) obtained by application of this Green's 
function to CDS pricing are summarised. 

In Section 5, similar expressions for the PV of other credit derivatives, specifically credit-contingent 
interest rate swaps (including with capped or floored Libor) and contingent CDS with an interest rate 
swap underlying are used in conjunction with those developed in Section 2 to assess the level of model 
risk associated with the uncertain parameter(s). Finally, in Section 6, we present some concluding 
remarks and a number of directions for possible future work. 


1 It should be mentioned that the author in his book eschews the idea that any one book should aim to be definitive. 
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2. Model Risk Methodology 

2.1. Previous Work 

A number of authors have previously visited the question of what constitutes an appropriate 
methodology for the quantification of model risk in pricing financial derivatives. In his pioneering 
work on the subject, Cont (2004) proposes two approaches. In the first, a family of plausible models is 
envisaged, each calibrated to all relevant market instruments and then used to price a given portfolio 
of exotic derivatives. The degree of variation in the prices that are observed provides a measure of the 
intrinsic uncertainty associated with modelling the price of the portfolio. A second approach, taking 
account of the fact that not all models are amenable to calibration to market instruments, compares the 
models by penalising them for the pricing error associated with calibration instruments. The pricing 
errors for multiple instruments can be combined using various choices of norm, giving rise to a number 
of possible measures of model risk. 

While intuitively attractive, neither of these approaches appears to have been adopted by 
practitioners. This is likely a consequence of the cost of implementing multiple models and re-pricing 
under them. Financial institutions usually only have very few models implemented, often just one, 
capable of pricing a given exotic option. Furthermore, regulatory pressure has recently been towards 
standardising pricing of financial derivatives by restricting or even reducing the size of the set of available 
models, which mitigates against the adoption of the kind of approach envisaged by Cont (2004). 

More recently, Glasserman and Xu (2014) have proposed an alternative approach based on 
maximising the model error subject to a constraint on the level of plausibility. The approach starts from 
a baseline model and finds the worst-case error that would be incurred through a deviation from the 
baseline model, given a precise constraint on the plausibility of the deviation. Using relative entropy 
to constrain model distance leads to an explicit characterization of worst-case model errors. In this 
way, they are able to calculate upper bounds on model error. They show how their approach can be 
applied to the problems of portfolio risk measurement, credit risk, delta hedging and counterparty risk 
measured through credit valuation adjustment (CVA). 

Although this approach has the attraction of a rigorous definition and, according to the authors, 
is amenable to convenient Monte Carlo implementation, it has the disadvantage that an entropy 
constraint specified a priori is not the sort of concept that risk managers are likely to be comfortable 
with in defining or expressing risk appetite. However, it is central to the whole approach. Furthermore, 
the approach has the disadvantage that it probably offers too much laxity in allowing the joint 
probability distribution function governing risk factors to vary freely subject only to the entropy 
constraint. Many of the perturbed distributions, including those giving rise to worst-case errors, would 
likely be deemed "unrealistic" by practitioners for reasons that cannot easily be encoded through 
entropy considerations. An approach that allows the user to be more specific about what is believed to 
be "known" and with what degree of certainty using a parametrisation more closely related to market 
variables would probably be preferred. 

For example, the consensus among practitioners might be that the "best" interest rate model 
would be somewhere between a normal and a lognormal process. However, under the proposal of 
Glasserman and Xu (2014), if a Hull-White (normal) model were chosen as the baseline, deviations 
towards lognormal and away from it would be penalised equally. However, we are really only 
interested in assessing the impact of the former. 

In his review of the subject Morini (2011), while lamenting the paucity of the model risk literature, 
comes down against excessive use of mathematical formalism and numerics that can serve to obscure 
the all-important link between specific modelling assumptions and the variability of prices that can 
emerge therefrom and advocates a middle path between that and "formal compliance or simple 
techniques to produce numbers that are acceptable to put in reports, but lacks [sic] the quantitative 
approach that would be needed to understand models deeply." A useful insight into the sorts 
of techniques that are currently used to put numbers into reports has recently been provided by 
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Joshi (2017). A more forward-looking perspective on the strategies that are being developed and 
implemented for model risk management in financial institutions is provided by Crespo et al. (2017). 

In the present work, we look to build to some extent on the basic philosophy of Cont (2004), 
but simplifying the methodology in compliance with the advocacy of Morini (2011) and avoiding the 
prohibitive cost of implementing multiple numerical models. We shall suggest that key to making 
progress is the ability to assess, at least to a good approximation, the impact of more advanced model 
features without necessarily having to implement them explicitly in a fully working model. 

To this end, we propose that asymptotic analysis that has, in the author's view, been under-used 
in risk management offers a fruitful way forward, certainly in the context of credit derivatives pricing 
with which we shall mainly be concerned here. To this end, we draw extensively on the asymptotic 
analysis of the impact of rates-credit correlation on CDS pricing in Turfus (2017a) and on contingent 
CDS in Turfus (2017b). A key advantage here is that the output of the model risk assessment process is 
in the form of analytic formulae rather than numerical routines. The relative transparency of the former 
will often furnish insight about which model configurations and inputs give rise to the greatest degree 
of uncertainty in terms of output prices. The alternative approach that is typically followed in model 
risk analysis is to sample the phase space of all possible market conditions, model parameterisations 
and product characteristics in a generally rather unsystematic way and infer conclusions from such 
data as are gleaned in the process. An obvious drawback of this approach is that there are no guarantees 
that the worst cases where the potential discrepancies are largest are uncovered; or how representative 
the numbers in the dataset are of the phase space that has been "sampled". 

We suggest that a further reason why the pricing model risk methodologies that have been 
proposed in the literature appear to have had limited traction in financial institutions lies in the 
structure of those institutions. Specifically, it is an issue that authority to make judgements in relation 
to which models are used, how they are configured and calibrated and how model risk is assessed 
is inevitably devolved across multiple functions: front office, market risk, model validation, finance, 
CVA desk, etc. For that reason, we propose that a simplified approach, whereby questions about how 
price-relevant model parameters are assigned and how the uncertainty associated with these is assessed, 
can be considered separately from issues around how their values and assumed uncertainty levels 
(or distributions) impact on prices. Our focus in this paper will by choice be on the latter. We note in 
this context that a further practical advantage of our proposed approach over the more mathematically 
sophisticated approach of Glasserman and Xu (2014) is that the intuition of practitioners tasked 
with the responsibility to assess parameter uncertainty levels is better attuned to more immediate 
parameters like correlation than the more abstract concept of entropy. 

2.2. Proposed Framework 

We formally state the problem we are looking to address as follows. Consider a model A4 (s, p ) 
which we wish to use as the basis for pricing a portfolio <E> containing derivatives D/,, k — 1, 2, ..., m. 
Here, s = (sj, S 2 , ..., s n ), with s, the value of a credit spread (suitably defined) associated with 
maturity T,, i — 1,2,... ,n, and p the correlation between rates and credit, the appropriate value of 
which is unknown and furthermore not readily ascertainable from market data, or else not in practice 
modelled. We wish to consider and indeed quantify, to a good approximation, the dependence of the 
portfolio price on the correlation parameter. 

To this end, we consider the calibration at time t = 0 of credit spreads in the model to a vector of 
market prices p — (pi, pi, ..., p n ) for calibration instruments {X\, X 9 , ..., I,,}. Let us express the 
result of such calibration for an assumed value of p formally as 

s = f{p;p) (1) 

for some / : [—1,1] —t M' ! . Since the market instrument prices are considered fixed for the time t — 0 
of interest, we shall for convenience generally omit the explicit dependence on p, in the following. 


Int. J. Financial Stud. 2018, 6, 39 


5 of 20 


in particular writing component-wise simply s, — fi{p). There may be other market instruments to 
which the model is calibrated, but, if the generated prices of these are not sensitive to p (or only very 
weakly so), we need not consider them explicitly in our analysis here. 

Let us then denote the price calculated for derivative D using model A4 (s, p), thus calibrated by 
V{D; Ai(s,p)). Let us further introduce the shorthand notation 

V k (p) — V(D k ;M(f(p),p)). (2) 


From the asymptotic analyses that will be described below, we see that the dependence of credit 
derivatives prices on p tends to be well captured as a linear function thereof. 2 An appropriate measure 
of the model risk associated with pricing the derivative portfolio with an assumed value p — po is on 
this basis obtained by use of the linear approximation: 


k=l 

with Ap an estimate of the level of uncertainty or inaccuracy associated with the representation of the 
parameter p. Performing the required differentiation on (2), we see 



dV k (p) 

dp 


P=P o 


dV(D k ;M(s;p)) 

dSj 


fi(p) + 


dV(D k ;M{s;p)) 

dp 


s=f(p),p=Po 


(4) 


We seek a convenient practical means of determining f- (po)- To that end, we note that the z’th 
calibration condition can be expressed as 


V(2pM(s;p)) = 


(5) 


leading by the same token to 


^dV(lp,M(s;p)) fU ^ , dV(I i; M(s;p)) 

^ - di- - h {po) + -^- 




dp 


— 0. 


s=f{p),p=Po 


( 6 ) 


Here, we assume that the model is bootstrapped by applying the calibration conditions in order of 
maturity, whence there will be no dependence of the price of X, on Sj for j > i. The partial derivatives 
in this equality can be computed conveniently by use of the asymptotic modelling approach described 
below. From this, the values of f- (po) can be inferred recursively. Substituting in (4) and then (3) 
gives rise to our representation of the model risk. On the basis of our assumption of approximate 
linear dependence of prices on p over the range of interest, we propose that our evaluations of /• (po), 
with typically p o = 0, can consistently be used in place of /■ ip). 

Our suggestion here is that, if we can derive analytic approximations to instrument prices taking 
into account the uncertain model parameters, this opens the way to obtaining analytic representations 
of the partial derivatives in (4) and so to obtaining an estimate of the model risk more conveniently 
and in a more transparent form than otherwise. We will illustrate our approach with some examples 
from the credit derivatives area, with which the author is most familiar. 


2 For other choices of model parameter than p, this will often still tend to be the case on the basis that the uncertain/unmodelled 
parameter will be of secondary importance; if this were not so, the consequent high degree of uncertainty introduced into 
pricing would compromise the utility of the pricing algorithm. 
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3. Two-Factor Asymptotic Model 

3.1. Underlying Processes 

Our modelling approach will be to represent the interest rate if and the credit default intensity Af 
(of a named debt issuer) as correlated mean-reverting short rate processes. In this respect, our approach 
is similar to that pioneered by Schonbucher (1999) who took both processes to be normal mean-reverting 
diffusions, in other words governed by the Gaussian short rate model of Hull and White (1990). 
Solutions were in his case found by constructing a two-dimensional tree. As was pointed out by 
Schonbucher (1999), it is a straightforward matter to extend his model to non-Gaussian processes. 

A number of authors have followed this suggestion taking the credit process to be lognormal, 
governed by a Black and Karasinski (1991) short rate model, which, although less tractable than 
a Gaussian model, ensures that credit spreads stay positive (and thus that survival probabilities are 
decreasing functions of time). Jobst and Zenios (2001) sought to price portfolios of bonds, modelling 
the credit spread for securities in a given rating class in this way, coupled with a Hull-White interest 
rate model, but also allowing rating class migrations to take place. A similar approach with only 
rates and credit default risk was used by Cortina (2007) to provide analytic solutions for the prices 
of defaultable bonds in the assumed absence of correlation, and by Pan and Singleton (2007) who 
considered the joint distribution of credit spreads and default loss rates implied by CDS market data. 

We will follow the latter authors in taking the interest rate process to be normal, as proposed 
by Hull and White (1990), and the credit intensity process to be lognormal, so ensuring positive 
intensities, following Black and Karasinski (1991). The correlation p between these two processes we 
take to be the uncertain model parameter of interest, although we could equally extend our framework 
to allow consideration of the credit mean reversion rate, or even its volatility as uncertain model 
parameters. We shall find it convenient to work with auxiliary variables Xt and i/f satisfying the 
following Ornstein-Uhlenbeck processes: 

dx t — —ot r x t dt + cr r (t) dW}, (7) 

dyt = -otxytdt + cr A (f) dWf, (8) 

where a. r/ a A € M + , <j r , cr^ : M + —> M + are piecewise continuous functions and dW}, dWf are correlated 
Brownian motions under the risk-neutral measure for t > 0 with 

corr (Wj 1 , W f 2 ) = p. (9) 

These auxiliary variables are related to the interest short rate jq and the credit default intensity Af, 
respectively, by 


r t — r(t) + r*(t) + Xt, (10) 

Af = (A(f) + A* (t))£(yt). (11) 

Note that here and below we use subscript-f to indicate stochastic processes; otherwise, all 
functions of f are assumed to be deterministic. Here, f(t) is the instantaneous forward rate 
(which we allow to take on negative values), A(f) the associated credit spread (see (15) below) and 
£ (Xt) exp { Xt — j [Xjf j is a stochastic exponential with [X]f the quadratic variation of a process 
Xf . Here, all functions and processes in (10) map [0, T m ] —>• M (although in practice r* (f) wil be strictly 
positive for t > 0), while those in (11) map [0, T m \ — » M + , where f = 0 is the "as of" date and T m the 
longest maturity date for which the model is calibrated. 

The interest rate model obtained in this way is of Hull-White type and the credit intensity 
model Black-Karasinski. It is well known that both of these models can be made risk-neutral by 
specification of suitable respective drifts, which determine in turn the configurable functions r* (■) and 
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A* (■) uniquely. The required form of r* (f) and A* (f) is in practice obtained by calibration of the model 
to ensure satisfaction of the no-arbitrage conditions set out below, assumed to apply for all te(o,T m }. 
As we shall see, it proves possible on this basis to determine these functions uniquely in asymptotic 
form under the conditions specified above. 


3.2. The No-Arbitrage Conditions 

The formal no-arbitrage constraints that determine the functions r* (f) and A* (f) are as follows: 

£ \e~ fo rsds 


Jo ( r s+A s ) ds 


= D(0,i), 
= B(0,f), 


under the martingale measure for 0 < t < T m , where 

„ \ — f / 2 ft s) ds 

D{t l ,t 2 ) =e J h w 

is the f]-forward price of the f 2 -maturity zero coupon bond and 

B( t 1 ,t 2 ) = e- I h {T{s)+I{s))ds 


( 12 ) 

(13) 

(14) 

(15) 


the corresponding risky bond price. We shall assume the bond prices can be ascertained at the initial 
time t — 0 from the market, whence we can view (14) and (15) as defining the forward rate r(f) and 
associated credit spread A(f), respectively. 

3.3. Derivation of Governing PDE 

We consider the general problem of pricing a cash security with maturity T whose payoff depends 
on xj. We will also look below at protection instruments whose payoff may depend on r and x T , 
where r is a stopping time in (0, T). We introduce the convenient shorthand notation that, for a process 
Xf and function / : M + —> M, 

£ x (f(t)X t ) := £(J(t)X t ))\ Xt=x , (16) 

in terms of which we can re-write (10) and (11) as r t = r(x t , t) and A t = A(y f , t), where 


r{x,t) r(t) + r*(t) +x t , 

\{y,t) := (A(t)+A*(0)£y(yt). 


(17) 

(18) 


Writing the price of the security at time t £ [0, T A t) as fj — f(xt,yt,t), we can infer by 
application of the converse of the Feynman-Kac theorem to (7) and (8) in the standard manner that the 
function f(x, y, t) satisfies the following backward diffusion equation: 


— + C-r{x,t)~ A (y, t) ) f(x, y, t) = 0, 


where 


- d 3 1 

L -0i r x- -a A y— + - 

dx J dy 2 


(19) 


af(l)A +2pcrr{l)aM ^ 




a 2 

dy 2 


( 20 ) 
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with in general lim fj = P(xj) for some payoff function P(x)? In the absence of closed form 

solutions to (19) and guided by the work of Hagan et al. (2005), Pagliarani and Pascucci (2011) and 
Horvath et al. (2017), we propose a perturbation expansion approach as follows. 

For both short rate models, we apply a 'low rates' assumption. To this end, we define 
small parameters 


0C r lm ^ 

fTm 

/ |r(f)l dt, 

1 0 

(21) 

i 

fTm _ 


• rp 

1 m 

•4-i 

<< 

__o 

(22) 


We assume that r(t ) and cr r (t) are 0(e r ), while A(f) is d(e^). The scaling of r*(f) and A*(f) is 
inferred as part of the calculation. We presage our conclusions by writing 

r*(0~7 io(t), 

A*(0~7i,i(0 +75,2(4 

with 7 *j(t) — 0{e r l e\i). We rewrite (19) as 

(Jj +C~r{t) - A(f) -tp e (x,y,t)J f(x,y,t ) = 0, (25) 

where 


(23) 

(24) 


^ e (x,y,f) := h(x,t) +g{y,t), 

(26) 

h(x, t) := r(x,t) -r(t), 

(27) 

g{y,t) := A(y, t) — A(f). 

(28) 


We take advantage of the assumed smallness of <p £ (.) to seek a Green's function solution C ( x, y, t ; £, ?/, v) 
for (25) as a joint power series in e r and e\, asymptotically valid in the limit as these two parameters 
tend to zero. The details of this calculation are provided in Turfus (2017a) with the conclusions 
summarised in Appendix A below. 

4. CDS Pricing 

We next consider how we can use our Green's function to price a credit default swap (CDS) 
analytically under an assumed rates-credit correlation, again following closely Turfus (2017a). 
Although this is a vanilla instrument, its use in calibration means that it is nonetheless important to 
have analytic formulae. 

4.1. Fixed Coupon Leg 

If, as proposed, the risky discount factors B (t \, t 2 ) are assumed known, a coupon payment made 
for a payment period f,] with coupon c and f, > 0 can be straightforwardly priced as 

PV Cou P on= cB (.°'ti)Ai, (29) 

with A, the relevant year fraction. 


We shall henceforth write for simplicity f(x, y, T) as a shorthand for lim f{x,y,t). 


3 
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4.2. Protection Leg 

Turfus (2017a) shows how the price of a protection leg can be derived by solving 
a nonhomogeneous version of (25) with the forcing function — (1 — R)A(ij, t) on the r.h.s., where 
R is the assumed recovery level of the referenced debt. The result obtained is, in the present notation. 


rT r oo rco 

PVp rot ~(l-R )/ / / A(ri,v)G{Q,Q,Q-,^,r],v)d^drjdv 

~ (1 — R) B(0,v) (A(v) + AA(v))dv (30) 

per unit notional with 0(e r (e r + ey) ) relative error, where 

AA{v) := - A(v) [ e~ ar ^ v ~ u) I rA (n) du, (31) 

with x(') given by (A28) and I rA (■ ) by (A9), provides an 0(e r e,\) adjustment to the leading order 
result. Here, the first term in the expression for AA(-) comes from the application of Go,o to A* (■) and 
the second from the application of Gj o to A(-). Note that, in the absence of correlation, AA(-) = 0 and 
the value of protection is as given under the assumption of deterministic rates. 

4.3. Calibration to CDS Market 

If we consider our model to be calibrated to risky bond prices, the calibration is at this stage 
completely specified, at least to second order accuracy. In particular, we see that f' (p) = 0 in (4), 
simplifying our task. 

Alternatively if, as is often the case, the calibration is to a term structure of CDS rates, we can take 
the market prices p, associated with maturities T; to be zero and the associated market instruments to 
be ATM CDS. Let us further suppose that the function A(f) can be taken as piecewise constant between 
the T/, given say by 

A(i)=Af, t g Tj\. (32) 

with Tg = 0. We can then take the s, introduced in Section 2.2 above to be given by 


s i T i = £A j (T j -T h i) 
;=i 


(33) 


Inference of the f- (p) is then straightforward, but our task is simplified if we are willing to 
consider only the leading order impact of calibration, whence we can neglect the 0(e A ) indirect impact 
of the A, through the (risky) discount factors in favour of their O (e A ) direct impact in the context of 
default-driven payoffs. A straightforward calculation gives rise to the conclusion that 


T t - f 0 T ‘ B(0,u)AA(u)du 
pTi ^ B(0, u)du 


(34) 


with expected 0(e A ) relative errors, which is consistent with our use of (30). 4 Equipped with this 
additional information, we are in a position to assess the model uncertainty associated with other 
derivative types priceable by our model. 


4 


The errors can in addition be expected to approximate to near zero since the calibration swaps are assumed to be at the 
money, whence the (risky) discounting affects both legs almost equally. 
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5. Calculating Correlation Risk 


5.1. Interest Rate Swap Extinguisher 

An interest rate swap extinguisher is an interest rate swap where the cash flows are contingent on 
survival of a named debt issuer. We have already considered credit-contingent fixed flows in Section 4.1 
above. We now look to price credit-contingent Libor flows. The payoff at time f, for a payment period 
tj] is given in the notation defined in Appendix A below by 

Payoff = X ti (x ti _ 1 ,t i _i)~ 1 - 1 

~ + x ti _ 1 B*(t i - tf-O) - 1, (35) 

with errors = 0(ej :). The calculation for the PV of this Libor flow contingent on no default was 
performed by Turfus (2017a). This was found to be given by 


PV, 


(0 


Libor 


B(0,ti) 


f 1-AlW 
y D(fj_i, fj) 



with 0(e^ + eV) relative error, where 


A J A (v)(p^\v)dv, 

(i) 

ya w) / 7 ( u 'V)I r \(u Av)du, 


7 (u,v) 


e ax ( v u \ u < v, 
e -x r {u-v)' u>v 


(36) 


(37) 

(38) 

(39) 


and I r \ (•) is given by (A9). We here use the binary operators A and V to represent min and max 
respectively. In conclusion, the fair price of a payer extinguisher will be 

PV Extinguisher = E (pv£L - Pv£ upon ) (40) 

1=1 

assuming the payments are synchronised. (The extension of the calculation if payments are not 
synchronised is trivial.) 

A comparison of analytic calculations based on (40) against the results of a finite difference solution 
of the underlying PDE is reproduced from Turfus (2017a) in Figure 1 to illustrate a typical parameter 
dependence structure and to indicate the level of accuracy that is furnished by our asymptotic method. 
The swap extinguisher paid quarterly Libor + 100 bp spread and received a quarterly 400 bp fixed 
coupon against a swap notional of 100. The credit default intensity was taken to be 770 bp, with a local 
vol of 60% and a mean reversion rate of 0.3. The 10-year swap rate was taken to be 80 bp with a short 
rate local volatility increasing from 20 bp to 70 bp and a mean reversion rate of 0.25. As can be seen 
from the graph, the use of our linear approximation approach to the model risk is a good one, with the 
discrepancy between the two modelling approaches in all cases less than 0.1 bp of notional. 
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Figure 1. Impact of correlation on PV for a 5 y maturity interest rate swap extinguisher. 


It is from here a straightforward matter of differentiation to quantify the model uncertainty 
associated with the parameter p. For the coupon flows, there is no such dependency to leading order. 
For the Libor flows, we have 


dPV, 


Libor 


d P 


5(0,h)AL(*) 

pD(f,_ 1 ,f ; ) 


(41) 


and, again ignoring the higher order indirect impact of the A j through the (risky) discount factors, 
we obtain 


dAi 


B( 0,ti) r(fi v b'-i)AT y (i) 

' (p y A '(v)dv. 


D(ti-l,t{) JTj —! 


(42) 


7 

From (4), we infer that, if the uncertainty associated with p is Ap, the model uncertainty associated 
with an interest rate swap extinguisher calibrated to risky bond prices is 


Uncertainty sa Ap 

and, if calibration is to CDS rates: 

Uncertainty « A p 


^dPV, 

i=1 


(0 

Libor 

dp 


(43) 


N / n apy(0 dPV® 


i =1 V=1 


3Ay 


(44) 


with fj (p) given by (34). Notice that, in the latter case, the impact of calibration adjustment is such as 
to reduce the overall uncertainty (for either a payer or a receiver swap), so ignoring it would be to take 
a conservative approach. 

5.2. Contingent CDS 

We consider next a contingent CDS on an interest rate swap with 10 years to maturity, paying 
semi-annual Libor + 40 bp and receiving a quarterly fixed coupon of 250 bp. The result below is from 
Turfus (2017b). We look to calculate the cost of providing protection against default of the counterparty, 
in other words the counterparty value adjustment (CVA) associated with the payer swap position. 
The value of protection up to some horizon T will be governed by the nonhomogeneous version of the 
governing PDE: 

+ £~r(t) -A(f) -<p e {x,y,t)^j f(x,y,t) = -A(y,t)P def (x,t), 


(45) 
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with Pdef{x, r) the protection payoff in the event of default at time x, subject to the final condition 
f(x,y,T) = 0. For the swap defined above, we can write 


p def( x ' T ) = max 


(1 - R) E (v^ (x, t) - cAiV^(*,r)) ,0 
i=l V 


(46) 


with 1? the counterparty recovery rate, where the Vp' > represent the PVs of the Libor flows and the 
V F ' those of the fixed coupon payments associated with the respective payment periods. Using a first 
order approximation based on (A18), we have 

Vp'(x,f) = X t{ (x, t) 

~ (1 — — t)), (47) 


with errors = 0(e $), where B* (•) is given by (A12). Likewise, we have to the same level of accuracy 



(x,f) 


X f -l(x,f) -X f i(*,f), 

.mgh-X'M- 


t < U-i, 

fj_l < f < tj. 


D(f,*,-_!) (1 - xB*(f;_i - f)) - D(f,h) (1 - *B*(f; - f)), 


D(f;_i, f) 


-l 


f) 

(*f—*i-l) 


- D(f,h) (1 - xB*(h - f)). 


t < k-1, 
h_i < f < f,-. 


(48) 


The solution of (45) is by standard application of the Green's function expansion (Al): only the 
leading order term Go,o( p ) is needed for our purposes. We conclude following Turfus (2017b) that, 
with (9(e^ + e?) relative error, the cost of protection purchased at f = 0 on a payer swap is given by 


N 

protection ~ (1 ~ R) £ (/f ~ cA //f ) 
t=l 


(49) 


where 

rWt\T _ 

fp := j A(v)B(0,v)D(v,w)N(—d\(£*(v),v))dv, (50) 

/f := (D(ti- 1 , ti)" 1 - l) P + B*(t, - *i_t) J^ T J(v)B(0,v)D(v,ti) 

(7(ti-i, v)Ir X (v A t,_i)N(-d 1 (f* (»), D )) + Ir (* A ) N,( - rf l (^0^’)) 

V vM») 

with /,(■) and I r ^ (•) given by (A7) and (A9), respectively, and 


di(x,v) 

X-I r x(v) 

(52) 

y/b(p) 


:= mf{x | P de f{x,v) > 0}, 

(53) 


where the latter expression need only be calculated to leading order, whence Xt t , can be replaced by x 
in (48). We further suppose formally that p def{') — G(e r ) as e r —» 0 so that the stochastic effects have 
a non-trivial impact, which is the situation of interest to us. 

A comparison of (49) against the results of a Monte Carlo simulation is reproduced from 
Turfus (2017b) in Figure 2 to illustrate a typical parameter dependence structure and to indicate 
the level of accuracy that is furnished by our asymptotic method. The credit intensity in this case is 
640 bp so not particularly "small"; the local volatility was taken to be 70% with a mean reversion rate 
of 0.3. The interest rate market was as in Section 5.1. The contract provided protection on the full value 
of the swap (assuming no recovery) for six years in return for semi-annual coupon payments of 400 bp. 


^ dv, (51) 
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The notional was again taken to be 100. As can be seen, the use of our linear approximation approach 
to the model risk remains good, with the discrepancy between the two modelling approaches unlikely 
to exceed a few basis points of notional. 
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Figure 2. PV dependence of interest rate swap protection on rates-credit correlation level. 

It is again a matter of straightforward differentiation to obtain an expression for the correlation 
risk associated with this modelling approach. To that end, we note that the relative impact through 
d i(-) in the second term in (51) will be (D(e\(e r + ex)), so we ignore this effect for our leading order 
calculation. We further note that, from (31), the relative impact of correlation through calibration is 
0(e r ex) and furthermore likely to cancel between legs, so we ignore this too. We obtain 


Uncertainty ss A p 


N rUAT 

fi/. 


A(z’)B(0, v)D(v, tj) 

- 1 - cA 


df rA (u At,—i) 


dp 


VW) 

+ B*(ti - t i _ 1 )j(t i _ 1 ,v)N(-d 1 ^*(v),v)))dv 


(54) 


It may be suggested that the computational effort required here could become burdensome if 
N were large. However, the greatest computational effort will be involved in computing £* (v) and 
the associated cumulative normal. The values of the latter can be tabulated in advance for a range of 
v £ [0, T] then interpolation used in the integration. Furthermore, for T = iy, we can re-express 



iAk r t: 

E J 

i=l Jt l -1 


and factor the integrand into the product of a u-dependent term and an /-dependent term, the latter of 
which can be taken outside the integral. This means we must integrate numerically from 0 to T only 
once, which is comparatively little effort. This approach was used to good effect by the author in the 
computations described in Turfus (2017b). 
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5.3. Capped Libor Floivs 

Finally, we consider the impact of capping a Libor flow such as was considered in Section 5.1 
above at some level K > 0, in this case focusing only on the leading order contribution. The calculation 
is similar to that presented in Turfus (2017c) where the Libor payment was assumed to be in a foreign 
currency, but the results below are new and presented here for the first time. The payoff at time t; for 
a payment period [f;_f ( ] is given in previously defined notation by 

Payoff = min{X f, (xf ; j, f z _i) —1 — l,KAj} 

~ min{D(f i _ 1 ,f ; )- 1 (l + x ti _ 1 B*(t i - f,^)) - l,KA t } (55) 


with errors = O(e^). As in the previous section, we assume the payoff to be 0(e r ) as e r —t 0 so that the 
stochastic effects have a non-trivial impact in that limit. Because of the appearance of x t _\ in the above 
expression, we must first compute the PV as of f;_i- We obtain by straightforward application of our 
leading order Green's function Coy: 


~ B(h-i,f ! )min{D(f ; _ 1 ,h) 1 (1 + x ti l B*(ti - f;_ a )) - 1,KA,}. 


(56) 


To proceed, we define the value 

*• := P(W;>CW1-1 (57, 

S*ft-t|-1)) 

as the (asymptotic) representation of the value of Xt t 1 at which the cap K is hit. Applying our (leading 
order) Green's function again to the payoff at f,_, to obtain the PV at t — 0, we obtain 

/ CO nOO 

/ G 0/0 (0,0,0; £, ?/, fj_! t ,_i ) dl, dij 

-00 7—00 

/ OO /•CO 

/ G 0 ,0 (0,0,0; f, 7/, f;_i) 

-OO 7x* 

(D(f ! -_ 1 ,f;)- 1 (l + P*(h - ii_i)) - 1 - XAf) (58) 

Carrying out the required integrations, we conclude 

PV Ca PP edU b or=mW) 

- B(0,ti)( (D(f i _ 1 ,t i )- 1 - l) N(d 1 (**,f / _ 1 )) +fCA / N(-d 1 (x*,h_ 1 )) 

- D(f;_ 1 ,f;) _1 B*(f I ' - ti—i)^/lr(ti_i)N'(-di(x*,t i _i))) (59) 

with 69(e f + e,\) relative error. In a similar vein, for a Libor flow floored at K, we have 

PV FlooredLibor ~ B(0, f f ) [k^N^X*, *,•_,)) + (D(f 1 _ 1 , f,)” 1 - l) N(-dr (**, f^)) 

+ D(f;_ 1; tfi-'B'iti- k_ x )^ (f ; _i )N'(—di(x*, t f _! ))). (60) 

On this occasion, none of the terms involving d,(-) should be neglected in calculating the model 
uncertainty since they constitute the leading order impact of correlation. They furthermore impact 
only one leg, not both, so there will be no cancellation between legs as in the previous case. Carrying 
out the necessary differentiation, we obtain, for both the capped and the floored cases: 
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Uncertainty r^i A P 


B(0 /ti ) dIrA ^- l)N ' { ^!r l)) ( i;)~ 1 - 1 - KAi 

Op y/l r {tj_ l) y 




■ ( 61 ) 


Again, we ignore the second order impact of correlation through calibration. We mention for 
completeness that, if the cap or floor is strongly in or out of the money, the uncertainty resulting 
from (61) will be small and possibly comparable with that of the plain Libor. We may in that case look 
to compute the latter separately using (44). 

6. Conclusions 

We have proposed a framework for the quantification of model risk in credit derivatives pricing 
in circumstances where the correlation between rates and credit is either uncertain in its value or 
not included in the calculation. We considered in particular the cases of (a) an interest rate swap 
extinguisher, (b) a contingent CDS on an interest rate swap underlying, and (c) an extinguisher with 
capped or floored Libor flows. We derived explicit analytic expressions for the model risk as a function 
of the degree of uncertainty associated with the correlation, under an asymptotic assumption of the 
interest rate and the credit default intensity being small and taking into account the potential impact 
of correlation on model calibration. We propose that the results obtained are accurate enough for risk 
management purposes. 

Although the cases considered here involve rather simple modelling considerations, we suggest 
that the approach advocated has much wider application, including for parameters other than 
rates-credit correlation in other multi-factor models. In particular, it is possible to look at asymptotic 
modelling involving also the price of a spot underlying such as an equity, an FX rate or an inflation rate. 
These quantities could further be assumed to jump in value contingent on default. Modelling then 
requires a three-dimensional diffusion process (possibly four, since two interest rates may appear, either 
or both of which may be assumed stochastic). Pricing of defaultable FX swaps, quanto CDS, contingent 
CDS on FX or equity options, convertible bonds and contingent convertible (CoCo) bonds can all 
be handled and analytic expressions for model uncertainty obtained in the manner specified above. 
For some examples of asymptotic formulae amenable to such analysis, see Turfus and Schubert (2017), 
Turfus (2017c), Turfus (2017d), Turfus (2017e) and Turfus (2018). In cases where a possible jump at 
default is assumed, the model uncertainty associated with uncertainty in the expected jump size can 
also easily be obtained as an analytic expression. 

Much work has also been done by other authors using perturbation approaches to obtain analytic 
approximations for prices of numerous option types under local and/or stochastic volatility modelling 
assumptions. See, for example, Pagliarani and Pascucci (2011), who considered equity option pricing 
under a local volatility assumption, obtaining a perturbation expansion for the relevant Green's 
function much as we did here, and using it to derive asymptotic expressions for option prices. Their 
approach was applied also to Asian option pricing in Foschi et al. (2013) and extended, with the 
use of some Fourier analysis, to incorporate Levy jumps in the dynamics of the spot underlying in 
Pagliarani and Pascucci (2013). A review of a number of other papers which have presented asymptotic 
pricing formulae in recent years has been given by Turfus and Schubert (2017). Our model uncertainty 
methodology is potentially applicable also to the results of such work. An interesting prospect for 
future work would be to combine asymptotic modelling of stochastic rates and local-stochastic volatility, 
as was done by Funahashi (2015), and look at the resultant model uncertainty in options pricing. 
Conflicts of Interest: The author declares no conflict of interest. 
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Appendix A. Green's Function 

From the analysis of Turfus (2017a), we infer that the Green's function solution of (25) can be 
expanded as 

CO 00 

G(x,y,t;£,rj,v) = (Al) 

i=0j=0 

with Gi j (■) = 0(e r , e i \i)- We will for the present purposes be interested only in terms up to second 
order, with consequent 0{e^ -\- e \) errors in G(-). We will in all cases be interested in "free-boundary" 
Green's function solutions which tend to zero as x,y —> ±oo. The leading order Green's function 
solution subject to these conditions is straightforwardly deduced. It is given by: 

Go,o(x,y,t;£,rj,v) = B(t,v)z^-N 2 (g - xe~‘ Xr(v ~ t \t] - ye~ a ^ v ~ tS> }R(t,v)), t < v, (A2) 


where N 2 (x,yiR(t,v)) is a bivariate Gaussian probability distribution function with mean 0 and 
covariance matrix 


R(t,v) 


hk {t,v) h(t,v) 


(A3) 


with 


Ir(h,h)'-= [ e llXr{tl U) a}(u)du, (A4) 

J h 

h(h, h) := / 2 e~ 1CL ^ u) o\{u) du, (A5) 

Jti 

Irk(h,h) ■= p [ e~ i0 ‘ r+0 ‘^ < ' t2 ~ u) o- r (n)a / \(u)du. (A6) 

J h 

For future notational convenience, we also define 


Ir(t) l r (0,t), (A7) 

h(t) := 4(0,f), (A8) 

Irk(t) := lrk(0,t). (A9) 

Following Turfus (2017a), we deduce at first order: 

Gi,o(x,y,t}£,ij,v) = - (xB* (v - t) + Go, 0 {x, y, t; g, t], v) (A10) 


and 

Go r i(x,y,t-,£,ri,v) = ~ J Mh) (£ y (e _flA(tl_t) i/f).A/(f,f 1 -1) G 0 jo(x,y,t;$,T],v) dt lr (All) 
where we have defined 


1 _ p-0LrT 

B*{ t):= -, (A12) 

Oiy 

:= J”e- a '( v - u h r {t,u)du, 

M tll t 2 G 0r0 (x r y r t- r ^,ij,v) : = G 0 , 0 (^x,y,t;^,i] - e~ a ^ v ~ t2 ' l I A (t 1 ,t 2 ),i^ ■ 


(A13) 

(A14) 
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Proceeding to second order, we obtain 




t 2 ) a 


-Xr(t2-t) dX y 

G 0/0 (x,y, t ; dt 2 dti 


+ It “ 7 2 ,o(b))rffiGo,o(x,y,t;^! 7 ,w), 


(A15) 


Go,2{x,y,t;Z,y,v) = J' A(f i){s y [e f) y f ) M t , h -l) ^ A(f 2 ) (f y (e aA(f2 At fl ,f 2 - l) 

G 0 ,o(x,y, f;£,y, v)dt 2 dt 1 

+ [\{h)S y (e-^-%) f l ’j{t 2 )Sy ( c -«A(fe-0 yt ) 

(exp ^ e -“A(f 2 -fi)j A (f,t 1 )^ - lj Mt,t 1 Mt 1 ,t 2 G Q j 0 (x,y,t-,g,n,v)dt 2 dt l 
- J { lo, 2 (h)£y (e““ A ( fl “ f) y t ) M t ,t 1 G 0fi {x / y,t;^ii,v)dt 1 (A16) 


and 


G U(I ,J,,«,,,,) - J' (»-«.-<> + £mw (4 (r-'A-V) X.„., -1) 

G 0i0 (x,y,t;£,ri,v) dt 2 dti 

+ /’!(,,> ( £y («-«(--)„) -1) /; 

Go,o(x,y,t m ,g,r],v) dt 2 dt 1 

+ ft It ^ tl)£ y { e ^ tl t] y^ e ^ tl h)l r\(t,h)Mt 1 ,t 2 GQ i Q(x,y,t;l,,ii,v)dt 2 dt 1 

+ f t Hh)£y (e~ akih ~ t} yt) B*(v-t 1 )I r ^{t,t 1 )Mt,t 1 G 0l o(x,y r t^,i],v)dt 1 

-J 7i,i(h)£y(e~° l ' i( - h ~ t ' > yt)Mt,t 1 G 0r0 (x r y,t;£,ri,v)dt 1 . (A17) 


Note that the order of integration between f ] and f 2 has been reversed compared to Turfus (2017a). 
It is a straightforward application of Fubini's theorem to derive the alternative expressions from 
the above. 

Finally, determination of the unknown 7 *■(■) functions is achieved by calibration of our model 
consistent with the no-arbitrage conditions (12) and (13). We must consider the consistent pricing in 
the former case of a risk-free cash flow, and in the latter case of a risky cash flow, as we now show. 

Appendix A.l. Calibration 

Appendix A. 1.1. Pricing of Risk-Free Cash Flow 

The calculation for a risk-free cash flow in our model is very similar to that performed by 
Horvath et al. (2017) and essentially corresponds to taking the distinguished limit as e\ —> 0 then 
e r —> 0. The same result is naturally obtained, namely that fj = X T (x, f), where, with the convention 
that Fij(x,t) — 0(e r l exi), 


X T (x,t) ~ D(t,T) (1 - F 1/0 (x,f) +F 2 ,o(x,t)) 


(A18) 
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with O(ef) errors, and our Green's function gives rise to 

Fijo(x,t) = xB*(T-t), (A19) 

- \x 2 B*{T - f) 2 + j\l*{t,v) - 7l, 0 (v))dv. (A20) 

Of interest to us here is the conclusion that, setting x = y — t = Oin (A18), satisfying (12) above 
to second order accuracy requires us to choose 

72 ,o(f) = -T*(0,f), (A21) 

which is 0(e 2 ), whence, on carrying out the required integration and using (A4), we can re-express 

FyoiX't) = \(x 2 — I r (t))B*(T — t) 2 — I*(T — t). (A22) 

The second term here is the convexity correction associated with the chosen money market 
numeraire, which term noticeably vanishes both at t — 0 and at t — T when the PV is 
known deterministically 

Appendix A. 1.2. Pricing of Risky Cash Flow 

We continue by writing the price at time f of a risky (zero recovery) cash flow at time T as 

fj — Y T (x t ,y h t), noting that, in this case, P(x) — 1 and /J = Y T (0,0,0) — B(0,T). We look to 

derive the general functional form of Y T (.) implied by our model, and in the process to determine the 
conditions on A* (f) necessary to satisfy (13) above. Applying our second order Green's function to 
this problem, we conclude 

Y T {x, y, t) ~ B(t,T) (l-F 1/0 (x,t) - F 0/1 (y,t) + F 2 , 0 (x,t) + F u (x,y,t) + F 0/2 {y,t)) (A23) 

with 0(ef + e \) error, where the F ;/ o( x , 0 are as defined above for i — 1,2 and we deduce, following 

Turfus (2017a), 

%( y,t) = J t A(p) (£ y (e^v-^yt) ~ l) dv, (A24) 

F o,i(y,t) = 2 f o,i(y '0 “ J f 7o,2( v )£y dv 

+ J t T \(v)£ y j V Mu)£ y (exp (w^ ’~ u) h(u)) - l) dudv, (A25) 

Fi,i{x,y,t) = xB*(T-t)F 0il (y,t) - ^ 7{ 1 {v)£ y (e ~“ A dv 
+ J t T \{v)£ y J V dudv 

+ J J \(u)£y j e~“A c - u ) j rA ( M ) dudv. (A26) 

Setting x — y — t — 0, we find that the no-arbitrage condition Y T (0,0,0) = B(0, T) is satisfied by 
the expression in (A23) to second order accuracy iff we choose 

To, 2 (0 = A(t) A(m) (exp (e" aA ( f "" ) t A (u)) - l) du, (A27) 

7 = J* (A(f) e - fl Ah-“) + \(u)e-^- u ^ I rX {u)du, (A28) 

which are (9(e 2 ) and 0(e r e x ), respectively. This completes the calibration of our model to second order. 
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